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0^ ■ Abstract. We study the joint effect of tidal torques and dynamical friction on the col- 
g i lapse of density peaks solving numerically the equations of motion of a shell of barionie 
i~ I . matter falling into the central regions of a cluster of galaxies. We calculate the evolution 

o ■ 

1 of the expansion parameter, a(t), of the perturbation using a coefficient of dynamical 
friction r) c i obtained from a clustered system and taking into account the gravitational 
2 ' interaction of the quadrupole moment of the system with the tidal field of the matter of 
^ | the neighboring proto-galaxies. We show that within high-density environments, such as 
^ . rich clusters of galaxies, tidal torques and dynamical friction slow down the collapse of 

' low-j/ peaks producing an observable variation of the parameter of expansion of the shell. 

b : 

5t i As a consequence a bias of dynamical nature arises because high-density peaks preferen- 
tially collapse to form halos within which visible objects eventually will condense. For a 
standard Cold Dark Matter model this dynamical bias can account for a substantial part 
of the total bias required by observations on cluster scales. 



1. Introduction. 

The origin and evolution of large scale structures is nowadays the outstanding problem in 
Cosmology. In a hierarchical 'bottom-up' scenario, high density collapsed peaks cluster 
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and merge to forme larger structures. The density fluctuation field is often assumed to 
be locally isotropic, the amplitudes are Gaussian distributed with uncorrelated phases 
(Peebles 1980). Under these assumptions the power spectrum for a continous random field 
is basically the squared amplitude of its Fourier modes. If we assume, for computational 
convenience, that the random field 5(r) is periodic in some large rectangular volume V, 
we can define the Fourier trasform to be: 

5(k) = i J (5(r)exp(ikr)d 3 r (1) 
According to Wiener-Khintchinc theorem: 

<|4| 2 >=-^/£(r)cxp(zkr)d 3 r (2) 
and 

C(r) = ^ / < |4 1 2 > ex P Hkr)d 3 r (3) 

i.e. the power spectrum, P(k), is the Fourier trasform of the autocorrelation function 
and vice versa. Also the density fluctuation field, 8(r), can be obtained from the Fourier 
trasform of the power spectrum, P(k). In other words, on average the characteristics of 
the density field peaks, e.g., their mass distribution, peculiar velocities, etc., are com- 
pletely determined by the spectrum through its moments (at least during the linear and 
early non-linear phases of the collapse (Bardeen et al. 1986)). Moreover the isotropy con- 
dition imposes that all physical quantities around density peaks is, on average, spherically 
symmetric. However, actual realizations of, e.g., the density field distributions around the 
density peaks which eventually will give birth to galaxies and clusters, depart from spher- 
ical simmetry and from the average density profile, producing important consequences 
on collapse dynamics and formation of protostructures (Hoffman & Shaham 1985; Ryden 
1988; Heavens & Peacock 1988; Kashlinsky 1986, 1987; Peebles 1990). A fundamental 
role in this context is played by the joint action of tidal torques (coupling shells of mat- 
ter which are accreted around a density peak and neighboring protostructures (Ryden 
1988)), and by dynamical friction (White 1976; Kashlinsky, 1986, 1987, Antonuccio & 
Colafrancesco 1995 (hereafter AC), Del Popolo & Gambera 1996). Some authors (see 
Barrow & Silk 1981, Szalay & Silk 1983, Peebles 1990) have proposed that non-radial 
motions would be expected within a developing proto-cluster, due to the tidal interaction 
of the irregular mass distribution around them, (typical of hierarchical clustering mod- 
els), with the neighboring proto-clusters. The kinetic energy of this non-radial motions 
opposes the collapse of the proto-cluster, enabling the same to reach statistical equilib- 
rium before the final collapse (the so called previrialization conjecture by Davis & Peebles 
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1977, Peebles 1990). Non-radial motions change the energetics of the collapse model by 
introducing another potential energy term. One expects that non-radial motions produce 
firstly a change in the turn around epoch, secondly a new functional form for S c , thirdly 
a change of the mass function calculable with the Press-Schechter (1974) formula and 
finally a modification of the two-point correlation function. Recently Colafrancesco, An- 
tonuccio & Del Popolo (1995, hereafter CAD) have shown that dynamical friction delays 
the collapse of \ow-v peaks inducing a bias of dynamical nature. Because of dynamical 
friction under-dense regions in clusters (the clusters outskirts) accrete less mass with 
respect to that accreted in absence of this dissipative effect and as a consequence over- 
dense regions are biased toward higher mass (Antonuccio & Colafrancesco 1995 and Del 
Popolo & Gambera, 1996). Dynamical friction and non-radial motions acts in a similar 
fashion: they delay the shell collapse consequently inducing a dynamical bias similar to 
that produced by dynamical friction but obviously of a larger value. This dynamical bias 
can be evaluated defining a selection function similar to that given in CAD and using 
Bardeen, Bond, Szalay and Kaiser (1986, hereafter BBKS) prescriptions. 

The plan of the paper is the following: in §2 we obtain the total specific angular mo- 
mentum acquired during expansion by a proto-cluster. In §3 we calculate the dynamical 
friction force for galaxies moving into the cluster, taking account of the clustering. In 
§4 we use the calculated specific angular momentum and the dynamical friction force to 
obtain the time of collapse of shells of matter around peaks of density having v c = 2, 3, 4 
and we compare the results with Gunn & Gott's (1972, hereafter GG) spherical collapse 
model. In §5 we derive a selection function for the peaks giving rise to proto-structures 
while in §6 we calculate some values for the bias parameter, using the selection function 
derived, on three relevant filtering scales. Finally in §7 we discuss the results obtained. 

2. Tidal torques and angular momentum. 

The explanation of galaxies spins gain through tidal torques was pioneered by Hoylc 
(1949). Peebles (1969) performed the first detailed calculation of the acquisition of angular 
momentum in the early stages of protogalactic evolution. More recent analytic compu- 
tations (White 1984, Hoffman 1986, Ryden 1988a) and numerical simulations (Barnes 
& Efstathiou 1987) have re-investigated the role of tidal torques in originating galaxies 
angular momentum. 

One way to study the variation of angular momentum with radius in a galaxy is that 
followed by Ryden (1988a). In this approach the protogalaxy is divided into a series of 
mass shells and the torque on each mass shell is computed separately. The density profile 
of each proto-structure is approximated by the superposition of a spherical profile, S(r), 
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and a random CDM distribution, e(r), which provides the quadrupole moment of the 
protogalaxy. As shown by Ryden (1988a) the net rms torque on a mass shell centered on 
the origin of internal radius r and thickness Sr is given by: 

(Irl 2 ) 1 / 2 = >/30(yC?) [{a 2m (r) 2 ){q 2 m(r) 2 ) 

-(a 2m (r)^ m (r)) 2 ] 1/2 (4) 

where qi m , the multipole moments of the shell and a; TO , the tidal moments, are given by: 

(12m (r) 2 ) = ^3 M 2 , | fc 2 dfcP (fc) j 2 (fcr) 2 (5) 

(a 2m (r) 2 ) = 1 dfcP (fc) , ?1 (fcr) 2 (6) 



r 



(a2 m {r)q* 2m {r)) = —p h M sh J kdkP (fc) ji (fcr) j 2 (fcr) (7) 

where M s /j is the mass of the shell, ji(r) and j 2 (r) are the spherical Bessel function of 
first and second order while the power spectrum P(k) is given by: 

P(k) = Ak- 1 [ln(l + 4.164fc)] 2 

(192.9 + 1340fc + 1.599 x 10 5 fc 2 + 1.78 x 10 5 fc 3 + 

3.995 x 10 6 fc 4 r 1/2 (8) 

(Ryden & Gunn 1987). The normalization constant A can be obtained, as usual, imposing 
that the mass variance at 8h~ 1 Mpc, as, is equal to unity. Filtering the spectrum on 
cluster scales, Rf — 3h~ 1 Mpc, we have obtained the rms torque, r(r), on a mass shell 
using Eq. (4) then we obtained the total specific angular momentum, h(r,v), acquired 
during expansion integrating the torque over time (Ryden 1988a Eq. 35): 

1 /3^ 2/3 



r o t 5/2 r (1-cosg) 3 Mf) 
M. h °° Jo (tf-smtf) 4 / 3 /iW-/ 2 W| 



the functions /i(^), f 2 ('&) are given by Ryden (1988a) (Eq. 31) while the mean over- 
density inside the shell, S(r), is given by Ryden (1988a): 

3 f°° 

5{r,v) = — daa 2 S(a) (10) 
r Jo 

In fig. 1 we show the variation of h(r, v) with the distance r for three values of the peak 
height v. The rms specific angular momentum, h(r,v), increases with distance r while 
peaks of greater v acquire less angular momentum via tidal torques. This is the angular 
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Fig. 1. The specific angular momentum, in units of Mq, Mpc and the Hubble time, t , for 
three values of the parameter v (y — 2 solid line, v = 3 dotted line, v — 4 dashed line) and for 
R f = 3/i _1 Mpc. 

momentum-density anticorrelation showed by Hoffman (1986). This effect arises because 
the angular momentum is proportional to the gain at turn around time, t m , which in 
turn is proportional to S(r, v)~i oc z/ -3 / 2 . 

3. Dynamical friction. 

In a hierarchical structure formation model, the large scale cosmic environment can 
be represented as a collisionless medium made of a hierarchy of density fluctuations 
whose mass, M, is given by the mass function N(M, z), where z is the redshift. In these 
models matter is concentrated in lumps, and the lumps into groups and so on. In such 
a material system, gravitational field can be decomposed into an average field, Fo(r), 
generated from the smoothed out distribution of mass, and a stochastic component, 
F s toch{f)i generated from the fluctuations in number of the field particles. The stochastic 
component of the gravitational field is specified assigning a probability density, W(F), 
(Chandrasckhar & von Neumann 1942). In an infinite homogeneous unclustered system 
VK(F) is given by Holtsmark distribution (Chandrasekhar & von Neumann 1942) while in 
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inhomogeneous and clustered systems T / F(F) is given by Kandrup (1980) and Antonuccio- 
Delogu & Barandela (1992) respectively. The stochastic force, F stoc h, in a self-gravitating 
system modifies the motion of particles as it is done by a frictional force. In fact a particle 
moving faster than its neighbours produces a deflection of their orbits in such a way that 
average density is greater in the direction opposite to that of traveling causing a slowing 
down in its motion. 

Following Chandrasekhar & von Neumann's (1942) method, the frictional force which is 
experienced by a body of mass M (galaxy) , moving through a homogeneous and isotropic 
distribution of lighter particles of mass m (substructure), having a velocity distribution 
n(v) is given by: 

= -47rG 2 M 2 n(w) J log Ap (11) 

where log A is the Coulomb logarithm, p the density of the field particles (substructure). 
A more general formula is that given by Kandrup(1980) in the hypothesis that there are 
no correlations among random force and their derivatives: 



fW(F)F 2 T(F)d 3 F 

F = -t]v = -± V 2^ v 12 

2 < v 2 > 



where r\ is the coefficient of dynamical friction, T(F) the average duration of a random 
force impulse, < v 2 > the characteristic speed of a field particle having a distance r ~ 
(^p^-) 1 / 2 from a test particle (galaxy). This formula is more general than Eq. (11) because 
the frictional force can be calculated also for inhomogeneous systems when W{F) is given. 
If the field particles are distributed homogeneously the dynamical friction force is given 
by: 

(Kandrup 1980), where m a and n a are respectively the average mass and density of the 
field particles. Using virial theorem we also have: 

i N 2 ' 3 (14) 



< v 2 > M tot 1 



Gm a n a /3 m n 1 / 3 ^, 
where M to t is the total mass of the system, R sys its radius and N is the total number of 
field particles. The dynamical friction force can be written as follows: 

11/2 

I ok i L.12.V- 1 \ 



F = _^ = _i^^ log{l , 2W v3 } _i_ (15) 



where N = ^-R 3 ys n a and n ac = n a x a 3 is the comoving number density of peaks of 
substructure of field particles. This last equation supposes that the field particles gener- 
ating the stochastic field are virialized. This is justified by the previrialization hypothesis 
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(Davis & Peebles 1977). 

To calculate the dynamical evolution of the galactic component of the cluster it is nec- 
essary to calculate the number and average mass of the field particles generating the 
stochastic field. 

The protocluster, before the ultimate collapse at z ~ 0.02, is made of substructure hav- 
ing masses ranging from 10 6 — 1O 9 M and from galaxies. We suppose that the stochastic 
gravitational field is generated from that portion of substructure having a central height 
v larger than a critical threshold v c . This latter quantity can be calculated (following AC) 
using the condition that the peak radius, r p k{v > u c ), is much less than the average peak 
separation n a (y > v^ 1 / 3 , where n a is given by the formula of BBKS for the upcrossing 
points: 

, ^ , exp(^/2) 7 o 2 

5fF^ exp( -^^ /18)1 (16) 

where 7, i?* are parameters related to moments of the power spectrum (BBKS Eq. 4.6A). 
The condition r p k(v > u c ) < 0.1n a (v > ^ c ) -1 ^ 3 ensures that the peaks of substructure 
are point like. Using the radius for a peak: 



r pk = V2R* 



1 



1/3 



(17) 



.(l + ^o)(7 3 + (0.9/^)) 3 / 2 

(AC), we obtain a value of v c = 1.3 and then we have n a (y > v c ) = 50.7Mpc~ 3 (7 = 0.4, 
R* = 50Kpc) and m a is given by: 

1 f°° 

m a = — — — - / m p k{v)N pk {v)dv = 1O 9 M (18) 

(in accordance with the result of AC), where m p k is given in Peacock & Heavens (1990) 
and N p k is the average number density of peak (BBKS Eq. 4.4). Clusters of galaxies 
are correlated systems whose autocorrelation function can be expressed, in the range 
10h~ 1 Mpc < r < 60/i _1 Mpc, in a power law form: 

£ cc - (^r (19) 

with 7 ~ 1.8 and a correlation length, r 0iC ~ 25h~ 1 Mpc (Bahcal & Soneira 1983; Postman 
et al. 1986). The analysis of fair samples of galaxies gives for the galaxy autocorrelation 
function the expression: 

z gg = C-^r (20) 
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in the range Q.lh^Mpc < r < 10h~ 1 Mpc (r . g ~ 5/i _1 Mpc, 7 = 1.77 ± 0.03 (Peebles 
1980, Davis & Peebles 1983)). The description of dynamical friction in these systems need 
to use a distribution of the stochastic forces, W(F), taking account of correlations. In this 
last case the coefficient of dynamical friction, 77, may be calculated using the equation: 

r] = J d 3 FW(F)F 2 T(F)/(2 < v 2 >) (21) 

and using Antonuccio & Atrio (1992) distribution: 

w[f) = ^f r dkksm ^ A f^ t 22 ) 

where Af is given in the quoted paper. The function Af is a linear integral function of 
the correlation function £(r). As shown in Del Popolo & Gambera (1997) the effect of 
clustering is that to increase the effects of dynamical friction. 

4. Shell collapse time. Tidal torques and dynamical friction acts in a similar fashion. 
As known one of the consequences of the angular momentum acquisition by a mass shell 
of a proto-cluster is the delay of the collapse of the proto-structure. As shown by Barrow 
& Silk (1981) and Szalay & Silk (1983) the gravitational interaction of the irregular mass 
distribution of proto-cluster with the neighbouring proto-structures gives rise to non- 
radial motions, within the protocluster, which are expected to slow the rate of growth 
of the density contrast and to delay or suppress collapse. According to Davis & Peebles 
(1977) the kinetic energy of the resulting non-radial motions at the epoch of maximum 
expansion increases so much to oppose the recollapse of the proto-structure. Numerical 
N-body simulations by Villumsen & Davis (1986) showed a tendency to reproduce this 
so called previrialization effect. In a more recent paper by Peebles (1990) the slowing of 
the growth of density fluctuations and the collapse suppression after the epoch of the 
maximum expansion were re-obtained using a numerical action method. In the central 
regions of a density peak (r < 0.5i?/) the velocity dispersion attain nearly the same value 
(Antonuccio & Colafrancesco 1997) while at larger radii (r > Rf) the radial component 
is lower than the tangential component. This means that motions in the outer regions 
are predominantly non-radial and in these regions the fate of the infalling material could 
be influenced by the amount of tangential velocity relative to the radial one. This can be 
shown writing the equation of motion of a spherically symmetric mass distribution with 
density n(r) (Peebles 1993): 

ft ft 71 f) 

_ n{Vr) + _ n(u 2) + (2(vl) - ( V j)) _ + n(r)-K) = (23) 

where (v r ) and (v$) are, respectively, the mean radial and tangential streaming velocity. 
Eq. (23) shows that high tangential velocity dispersion ((u^) > 2(v 2 }) may alter the infall 
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Fig. 2. The time evolution of the expansion parameter. The solid line is a(t) for GG model; the 
dashed line is a(t) taking account only dynamicalfriction; the dotted line is a(t) taking account 
of the cumulative effect of non-radial motions and dynamical friction. 

pattern. The expected delay in the collapse of a perturbation, due to non-radial motions 
and dynamical friction, may be calculated solving the equation for the radial acceleration 
(Kashlinsky 1986, 1987; AC; Peebles 1993): 

dv r L 2 (r, y) , . dr ,„ 

^ = w-^-^ (24) 

where L(r, v) is the angular momentum and g(r) the acceleration. Writing the proper 
radius of a shell in terms of the expansion parameter, a(r i7 1): 

r{n,t) = na{ri,t) (25) 
remembering that 

M=—p b { n ,t)a 3 {r l ,t)rf (26) 

and that pb — |^ , where Hq is the Hubble constant and assuming that no shell crossing 
occurs so that the total mass inside each shell remains constant, that is: 
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Fig. 3. The time of collapse of a shell of matter in units of the age of the universe t for v — 2 
(dotted line) compared with Gunn & Gott's model (solid line). 



p(n,t) 



Pi(n,t) 



a 3 (r u t) 

the Eq. (24) may be written as: 



(27) 



+ 



4G 2 L 2 



da 



(28) 



(Pa _ H 2 {l + 5) 

~d¥ ~ 2a~ 2 ' H 4 (l + 6) 2 rj°a 3 '' dt 
The equation (28) may be solved using the initial conditions: (^|) = 0, a = a max ~ 1/5 
and using the function h(r, v) = L(r, u)/M s h found in §2 to obtain a(t) and the time of 
collapse, T c (r, v). 

In Fig. 2 we show the effects of non-radial motions and dynamical friction separately. As 
displayed non-radial motions have a larger effect on the collapse delay with respect to 
dynamical friction. In Figs. 3 -f- 5 we compare the results for the time of collapse, T c , for 
v = 2, 3, 4 with the time of collapse of the classical GG spherical model: 



(29) 



As shown the presence of non-radial motions produces an increase in the time of collapse 
of a spherical shell. The collapse delay is larger for low value of v and becomes negligible 
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Fig. 4. The time of collapse of a shell of matter in units of the age of the universe t for v — 3 
(dotted line) compared with Gunn & Gott's model (solid line). 

for v > 3. This result is in agreement with the angular momentum-density anticorrclation 
effect: density peaks having low value 1 of v acquire a larger angular momentum than high 
v peaks and consequently the collapse is more delayed with respect to high v peaks. 
Given T c (r, v) we also calculated the total mass gravitationally bound to the final non- 
linear configuration. There are at least two criteria to establish the bound region to a 
perturbation 5(r): a statistical one (Ryden 1988b), and a dynamical one (Hoffman & 
Shaham 1985). The dynamical criterion, that we have used, supposes that the binding 
radius is given by the condition that a mass shell collapse in a time, T c , smaller than the 
age of the universe io : 

T c (r,u)<t (30) 

We calculated the time of collapse of GG spherical model, T c q(t, v), using the density 
profiles given in Ryden & Gunn (1987) for 1.7 < v < 4 and then repeated the calculation 
taking into account non-radial motions obtaining T c (r, v). Then we calculated the binding 
radius, rb {v), for a GG model solving T c0 (r, v) < t for r and for several value of v, while 
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Fig. 5. The time of collapse of a shell of matter in units of the age of the universe t for v = 4 
(dotted line) compared with Gunn & Gott's model (solid line). 

we calculated the binding radius of the model that takes into account non-radial motions, 
rj,(i/), repeating the calculation, this time with T c (r, v) < to. We found a relation between 
v and the mass of the cluster using the equation: M = ^-r\ph- 

In fig. 6 we compare the peak mass obtained from GG model, using Hoffman & Sha- 
ham's (1985) criterion, with that obtained from the model taking into account non-radial 
motions. As shown for high values of v [y > 3) the two models give the same result for 
the mass while for v < 3 the effect of non-radial motions produces less bound mass with 
respect to GG model, decreases the effect of non-radial motions produces a decrease in 
the bound mass. 

5. Tidal field and the selection function 

Following BBKS we define a selection function t(y/v t ) which gives the probability 
that a density peak forms an object, while the threshold level, v tl is defined so that 
the probability that a peak form an object is 1/2 when v — v t . The selection function 
introduced by BBKS (Eq. 4.13), is an empirical one 

and depends on two parameters: the threshold v t and the shape parameter q: 
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Fig. 6. The mass accreted by a collapsed perturbation, in units of 1O 1j M0, taking into account 
non-radial motions and dynamical friction effect (dotted line) compared to Gunn & Gott's mass 
(solid line). 



t{v/v t ) 



(31) 



1 + [vjvtY 

If q — ► oo this selection function is a Heaviside function $(1/ — v t ) so that peaks with 
v > v t have a probability equal to 100% to form objects while peaks with v < v t do not 
form objects. If q has a finite value sub-^ t peaks are selected with non-zero probability. 
Using the given selection function the cumulative number density of peaks higher than v 
is given, according to BBKS, by: 



/•OO 

n p k = / t(v/v t )N pk (v)dv 

J V 



(32) 



where N p k{v) is the comoving peak density (see BBKS Eq. 4.3). A form of the selection 
function, physically motivated, can be obtained following the argument given in CAD. 
In this last paper the selection function is defined as: 



/•OO 



dS 



(33) 



where the function 
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p\S,(8)(r)] 



1 



exp 




\S-(S)(r)\ 2 



) 



(34) 



gives the probability that the peak overdensity is different from the average, in a Gaussian 
density field. The selection function depends on v through the dependence of 8(r) from 
v. As displayed the integrand is evaluated at a radius rut which is the typical radius of 
the object we are selecting. Moreover the selection function t{v) depends on the critical 
overdensity threshold for the collapse, 8 C , which is not constant as in a spherical model 
(due to the presence, in our analysis, of non-radial motions and dynamical friction that 
delay the collapse of the proto-cluster) but it depends on v. An analityc determination 
of 8 c {v) can be obtained following a technique similar to that used by Bartlett & Silk 
(1993). Using Eq. (28) it is possible to obtain the value of the expansion parameter of 
the turn around epoch, a max , which is characterized by the condition = 0. Using the 
relation between v and Si, in linear theory (Peebles 1980), we find 



where 8 co = 1.68 is the critical threshold for GG model and A G and /i(8) are given 
in Colafrancesco, Antonuccio & Del Popolo (1995) (Eq. 5, 6). In Fig. 7 we show the 
overdensity threshold in function of v. As shown, 8 c {v) decreases with increasing v. When 
v > 3 the threshold assume the typical value of the spherical model. This means, according 
to the cooperative galaxy formation theory, (Bower et al. 1993) that structures form more 
easily if there are other structures nearby, i.e. the threshold level is a decreasing function 
of the mean mass density. Known 8 c (v) and chosen a spectrum, the selection function 
is immediately obtainable through Eq. (33) and Eq. (34). The result of the calculation, 
plotted in fig. 8, for two values of the filtering radius, (Rf — 2, 3 h~ 1 Mpc), shows that 
the selection function, as expected, differs from an Heaviside function (sharp threshold). 
The value of v at which the selection function reaches the value 1 (t(v) ~ 1) increases 
for increasing values of the filtering radius, Rf. This is due to the smoothing effect of the 
filtering process. The effect of non-radial motions is, firstly, that of shifting t(v) towards 
higher values of v, and, secondly, that of making it steeper. The selection function is 
also different from that used by BBKS (tab. 3a). Finally it is interesting to note that the 
selection function defined by Eq. (33) and Eq. (34) is totally general, it does not depend on 
the presence or absence of non-radial motions. The latter influence the selection function 
form through the changement of 8 C induced by non-radial motions itself. 



6. The bias coefficient One way of defining the bias coefficient of a class of objects 
is that given by (BBKS): 




(35) 
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b(R f ) = + 1 (36) 



where (v) is: 

7 6» 



< 



Jo 



t{-)N pk {v)dv (37) 
ft 



1 — 7 : 

from Eq. (37) it is clear that the bias parameter can be calculated once a spectrum, P(k), 
is fixed. The bias parameter depends on the shape and normalization of the power spec- 
trum. A larger value is obtained for spectra with more power on large scale (Kauffmann 
et al. 1996). In this calculation we continue to use the standard CDM spectrum (f2o = 1, 
h = 0.5) normalized imposing that the rms density fluctuations in a sphere of radius 
%h~ x M-pc is the same as that observed in galaxy counts, i.e. a% — cr(Sh^ 1 Mpc) = 1. The 
calculations have been performed for three different values of the filtering radius (Rf = 2, 
3, 4 h~ 1 Mpc). The result of the calculation is plotted in table 1. As shown, the value 
of the bias parameter tends to increase with Rf due the filter effect of t(v). As shown 
t(u) acts as a filter, increasing the filtering radius, Rf, the value of v at which t(y) ~ 1 
increases . In other words when Rf increases t(v) selects density peaks of larger height. 
The reason of this behavior must be searched in the smoothing effect that the increasing 
of the filtering radius produces on density peaks. When Rf is increased the density field 
is smoothed and t(v) has to shift towards higher value of v in order to select a class of 
object of fixed mass, M. 



Bias 



R f {hr x M V c) 


b 


2 


1.6 


3 


1.93 


4 


2.25 



Table 1. Values of the coefficient of bias for different values of the filtering radius. 



7. Conclusions 

In this paper we have studied the role of non-radial motions and dynamical friction 
on the collapse of density peaks solving numerically the equations of motion of a shell of 
barionic matter falling into the central regions of a cluster of galaxies. We have shown 
that non-radial motions and dynamical friction produce a delay in the collapse of density 
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Fig. 7. The critical threshold, 8 c (v) versus v 

peaks having low value of v while the collapse of density peaks having v > 3 is not 
influenced. A first consequence of this effect is a reduction of the mass bound to collapsed 
perturbations and a raising of the critical threshold, S c , which now is larger than that 
of the top-hat spherical model and depends on v. This means that shells of matter of 
low density have to be subjected to a larger gravitational potential, with respect to the 
homogeneous GG model, in order to collapse. The delay in the proto-structures collapse 
gives rise to a dynamical bias similar to that described in CAD whose bias parameter 
may be obtained once a proper selection function is defined. The selection function found 
is not a pure Heaviside function and is different from that used by BBKS to study the 
statistical properties of clusters of galaxies. Its shape depends on the effect of non-radial 
motions and dynamical friction through its dependence on 5 c (v). The function t{y) selects 
density peaks higher and higher with increasing value of Rf due to the smoothing effect 
of the density field produced by the filtering procedure. Using this selection function and 
BBKS prescriptions we have calculated the coefficient of bias, b. On clusters scales for 
Rf = Ah^Mpc we found a value of b = 2.25 comparable with that obtained from the 
mean mass-to-light ratio of clusters, APM survey, or from N-body simulations combined 
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Fig. 8. The selection function, t(v), for Rf = 3h 1 Mpc (S c = 1.68, solid line; S c function of v, 
dotted line) and for 4h~ 1 Mpc (8 C = 1.68, short dashed line; S c function of v, long dashed line). 

with hydrodynamical models (Frenk et al. 1990). Besides, the value of the coefficient of 
biasing b that we have calculated is comparable with the values of b given by Kauffmann 
et al. 1996. This means that non-radial motions and dynamical friction play a significant 
role in determining the bias level. 
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